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ABSTRACT 

We consider a system of evolution equations for quark and gluon structure func- 
tions satisfying the leading-logarithmic behaviour due to both QCD collinear (LLQ 2 ) and 
infrared (LLl/x) singularities. We show that these equations leave undetermined an arbi- 
trary regular function of j in the Mellin-transformed weights. We consider the constraints 
resulting from energy-momentum conservation and from the decoupling of quark loops in 
the leading j-plane singularity. These constraints can be fulfilled without influencing the 
leading-log terms. As a particular consequence of the second constraint, the location of the 
leading singularity is determined in terms of the (LLl/x) and (LLQ 2 ) kernels. It leads 
to a value significantly lower than the LLl/x evaluation, while remaining at j > 1, and 
compatible with the behaviour of structure functions observed at HERA. 
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The recent results on quark and gluon structure functions at Hera have paved the way 
for a reconsideration of QCD predictions. The observed behaviour of the proton structure 
function F2 (x, Q 2 ) and some indications on the gluon structure function Fq (x, Q 2 ) in the 
range {lO" 4 < x < 10" 2 , 8 < Q 2 < 60 GeV 2 } are characterized W by a rapid rise at small 
x which is qualitatively compatible with the predictions of the resummation of leading 
log(l/x) contributions (hereafter denoted LLl/x) of the perturbative expansion, i.e. with 
the Lipatov singularity ^ (BFKL). There exists quantitative studies of structure functions 
including the BFKL singularity ! 3 L 

On the other hand, new tests of the celebrated Altarelli-Parisi evolution equations^ 
(DGLAP) are now possible in a much larger Q 2 range . These equations correspond to the 
resummation of the leading log(Q 2 ) terms (denoted LLQ 2 ) of the perturbative QCD ex- 
pansion. Thus, HERA represents a unique apparatus for testing the QCD theoretical tools 
in kinematical domains where the perturbation expansion has to be resummed. As well- 
known, this is due to the appearence of collinear and infrared (in the infinite-momentum 
frame) singularities in the perturbative theory, implying large logarithms in the effective 
coupling constant. Hopefully, a precise comparison between experiment and theory will 
help understanding the yet unknown non-perturbative regime. In this context, a unified 
description of the LLl/x and LLQ 2 evolution equations for quark and gluon structure 
functions is highly desirable. 

Two different approaches have already been proposed, each one with its own advan- 
tages and inconveniences. First, it was remarked! 5 ' that Feynman diagrams for multi-gluon 
emission which contribute to Fq (x, Q 2 ) are characterized by a common angular ordering 
property in the whole x-range, leading after resummation to the BFKL singularity at small 
x and to the DGLAP equations elsewhere. This allows one to write a unique equation for 
the whole x ranged. This approach leads to fruitful Monte Carlo simulations ^ based on 
the multi-gluon diagrams but an explicit solution of the evolution equation! 6 ] itself has not 
yet been found. 
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In a second existing approach, explicit solutions have been proposed in a different 
context t 8 '. Using a constructive procedure starting directly from the structure functions, 
a system of evolution equations can be written following the scheme: 

(LLl/x) + (LLQ 2 ) - (DLL), (1) 

which indicates that the sum of LLl/x and LLQ 2 contributions is corrected for double- 
counting by substraction of the double-leading-logarithmic terms (DLL). This led the 
authors^ 8 ' to propose a set of explicit equations for the quark-gluon system. However, as 
we shall see below, we find that the realization of the scheme (1) in terms of an explicit set 
of evolution equations leads to a solution which differs from that proposed in Ref. [8] . On 
the one hand, the set of equations obtained there leaves undetermined a regular function 
related to higher order contributions. On the other hand, the realization of (1) has to 
obey specific constraints which should be satisfied by the resulting equations, namely i) 
energy-momentum conservation, and ii) the decoupling of quark loops in the BFKL kernel. 
These constraints are not satisfied by the system of Ref. [8] . 

The goal of this paper is to provide an explicit realization of the formal scheme (1), 
obeying the constraints. We thus provide an explicit solution of the fundamental equation 
written in Ref. [6], while extending it to the quark-gluon system. 

The generic set of equations we propose following the scheme (1) reads: 



dlnQ^ =® fG{3,Q)-QofG{3,Qo)+J Q , —^ d ® 

[{{u G + *) (j) Q 2 K (Q 2 , Q' 2 ) - V(j)}f G {j, Q' 2 ) + ^(j)f s (j, Q' 2 )} 

dlnQ^ =® fs{3,Q)-Qofs{3,Qo)+J Q , ~^^ d Q 

[MJ)fs (j, Q' 2 ) + 2n F ^(j)f G (j, Q' 2 )] , (2) 

where the inhomogeneous terms of the system of equations is implied by the compatibility 
with the DGLAP evolution equations near the threshold Qq. ^(j) is an arbitrary function 
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provided it remains regular at j = 1, and where 

a s {Q 2 ) = l/(MnQ 2 /A 2 ); {v G , <f>% v F , 0g} (j), 

are respectively the first-order QCD coupling constant, and the set of Altarelli-Parisi 
weights written with the usual conventions ^ . Here K [Q 2 , Q' 2 ) is the BFKL kernel ^ . Note 
that the specific equations of Ref.[7] are recovered with the choice \1/ = 4N c /(j — 1) — vq- 
As we shall now demonstrate, this choice is not required by (1) and moreover it contradicts 
energy-momentum conservation and the decoupling of quark loops in the BFKL kernel. 

Let us briefly show how the system (2) verifies both LLl/x and LLQ 2 behaviour, 
while leaving undetermined the regular function ^(j). It is well known^ that the kernel 
K (Q 2 ,Q' 2 ) reduces simply to l/Q 2 when is retained only its contribution to the LLQ 2 
terms. In that precise case, it is easy to realize by substitution in (2) that the function 
^(j) disappears from the system. Integrating by part Eqns.(2), one recovers the ordinary 
DGLAP evolution equations up to negligible terms. 

The LLl/x reduction of the system (2) is recovered by selecting the most singular 
contributions in the complex j-plane. This is obtained first by considering the equation 
with only the BFKL kernel K (Q 2 ,Q' 2 ) . Then one has also to select the most singular 
contribution from the prefactor ug + ^ ~ 4AT C / — By this way, one finds back the usual 
BFKL kernel. Again, the function ^(j) remains undetermined, provided that it is chosen 
regular at j = 1. By this discussion, we have proved that the requirement of equation (1) 
has a large class of solutions parametrized by the regular but arbitrary function 

Now, let us introduce the constraints on the system of equations (2) which cannot 
be satisfied by the leading-log terms only . The first one is due to energy-momentum 
conservation which reads: 

VQ 2 , {F G (j, Q 2 ) + F S (j, Q 2 )) \ j=a = 1 => f G (2, Q 2 ) + f s (2, Q 2 ) = 0. (3) 

Inserting the constraint (3) (for both Q 2 and Qq) in the system (2) taken at j = 2, one 
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gets: 

nQ " a s {Q' 2 ) dQ' 2 



= J q2 ^ [Q J Q (*(2) + i/ G (2)) (K (Q 2 , Q' 2 ) - 1/Q 2 ) , (4) 

where one uses the values of the Altarelli-Parisi kernels at j = 2, namely 

f 2 8 8 1 1 

Then, if the kernel if (Q 2 , Q /2 ) reduces to l/Q 2 the identity (4) is trivially fulfilled. Since 
it is not the case beyond the double-leading-log approximation, energy-momentum con- 
servation requires the relation \I>(2) = §np. This is the first condition on the previously 
undetermined function 

Let us now consider solutions of the system (2). In the present paper, we shall restrict 
ourselves to the conventional case of a fixed coupling constant as = a. Indeed, this is 
the only case for which the BFKL singularity has been derived from the diagrammatic 
expansion! 2 ' . Following a common procedure, usual in the LLl/x case! 2 ' 8 ' 10 ], we write the 
equations for fa,s m terms of their inverse Mellin transforms, namely: 

f G ,s(j,Q 2 )= j c ^^s{jn){Q 2 y-\ (5) 

where C is a contour in the complex-7 plane leaving all singularities of the integrand at its 
left. Inserting the Mellin transform (5) in the system (2), and after some manipulations, 
one finds the following matrix equation: 



'l-*{(^ + *)w(7)-f} 



4-7T 7 



a 2nF^ G ^ a_ vp 

4ir 7 4-7T 7 

dlnT 





f Vg) 


H 


f ¥>0G \ 




\<ps) 







(6) 



where a; (7) = 2-0(1) — ipd) — ^(1 — 7)5 with ip^) = ^ , is the well-known eigenvalue of 
the BFKL kernel 1 . Note that the vector of the right-hand side (6) includes the boundary- 



1 The scale dependence of the system has been merged into the right-hand side of eq.(6), 
since it is not relevant for the computation of the BFKL singularity at fixed coupling. 
However a treatment with a running coupling constant would imply to take into account 
the threshold Q 2 ,, see for instance Ref.[10]. 
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value conditions on the system (2). 
j-plane singularity. 

By matrix inversion one gets: 



It is usually assumed not to contain the dominant 




1 



\ £ 2n F ^ G 1 - A { W (7) (i/ G + *) - f } y \<p 0F 





f V?0G \ 







(7) 



where 



^(J',7) = 1 - £ |^(7) ("G + *) + ^ - *) [> , 



neglecting contributions of order a 2 . 

In order to obtain the leading behaviour of structure functions at low x, we perform 
a double inverse-Mellin transform on the system (7), first in the variable j, then in the 
variable 7 : 

fa,s Q 2 ) - / ^~ 1)lnQ2 / ^ e °" 1)lnl/X <PgM 7)- (8) 

The first integral of Eq. (8) (in the complex j-plane) is given by the residue of the pole at 
j = jp{l) where T>(j p , 7) = 0. Then the second Mellin transform (in the complex 7-plane) 
is performed by a saddle-point method. In the domain of validity of the BFKL singularity, 
one assumes <C 1, and the saddle-point is at 7 = 7 C corresponding to dj '^ c ' > = 0. 

One finds: 

4n 1 

(pole) — = (u G (j p ) + * (j p )) w( 7 ) + (u F {j p ) - * (j p )) - 
a 7 

(saddle point) = (y G (j) + * (j)) J ( 7c ) - \ (v F (j) - (])) , (9) 

7c 

where j = j p {^ c )- The resulting structure functions behave like: 

/cs^Q^^IQ 2 ] 76 -^-^ 1 ), (10) 

Now, we are able to express the second constraint on the system of equations (2), 
namely the decoupling of quark loops from the BFKL singularity kernel. Let us first 
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require that the dominant behaviour (10) agrees with the theoretical determination of 
the BFKL singularity ^ 2,8 \ It corresponds to the well-known "critical" value of 7, namely 
7 C = 1/2, with u/ (7 C ) = and uj{^ c ) = 41n2. From inspection of the system (9), this 
implies the following relations: 



* G) = Q) 



VG G) + V F G) 



Act log 2 
An 



-1 



(11) 



We note that the first of relations (11) implies the cancellation of the quark contribu- 
tion vp — ^ in the pole location, see (7). Indeed, it is known! 2 ! that a loop contribution 
for a particule of spin s is of order x~^ 2s ~ 1 ^ at small x. Thus, quark loops are negligeable 
in the BFKL kernel. 

The results of Eqs.(ll) have some interesting consequences on the location of the 
BFKL singularity in the j-plane. 

If we retain only the dominant behaviour vq ~ 4^-, when j 1— > 1, we recover the 
usual determination of the BFKL singularity, in the absence of LLQ 2 contributions. The 
coupling to the LLQ 2 terms, which is allowed by our unified description, is responsible 
for the modified formula (11). As an illustration, let us for instance fix the value of a 
such that A In 2 ~ .5. This value is often considered as fixing the "bare" BFKL 
singular ityt 3 '. Using this value and fixing N c = uf = 3, one gets from the second equa- 
tion (11): j — 1 = 0,313. Quite interestingly, this value is in better agreement with those 
phenomenologically obtained from Hera experimental results, if analyzed directly in terms 
of a BFKL singularity^ 1 '. Note that, with these input numbers, one gets for the function 
* the two constraints #(J) = v F (j = 1.313) ~ -1.44; *(2) = u G (2) = 2. The kine- 
matical constraint of energy-momentum conservation decreases the intercept of the BFKL 
singularity^ 11 !, through non-leading terms in the LLl/x expansion. However, this con- 
straint at j = 2 does not allow a determination of the modified location of the singularity. 
In our derivation, see equation (11), we obtain a precise theoretical prediction for this 
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shift, based on a different constraint, e.g. the decoupling of the quark loops from the 
BFKL kernel. 

Note that in all realistic cases (values of tif, flavour thresholds, etc.), vq {j)+vf (j) < 
AN/ (j — 1), for j > 1. Thus the decoupling constraint always decreases the BFKL intercept 
J-l- 

In conclusion, we have proposed QCD evolution equations combining the leading- 
logarithmic contributions in both the x and Q 2 variables and satisfying constraints, namely 
energy-momentum conservation and quark decoupling from the BFKL singularity. Our 
results extend the QCD predictions beyond the leading-logarithm contributions. In par- 
ticular, we find the location of the dominant singularity of structure functions at small-x 
and fixed coupling constant, which is shifted down with respect to the LLl/x prediction. 
Other consequences of this improved leading-logarithm scheme should be explored. In 
particular, it would be interesting to see in which way the system (2), with an adequate 
function ^(j), gives an explicit realization of the gluon resummation of Ref. [5]. 
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